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Abstract 

We investigate the quintom model of dark energy in the generaUzed case where the cor- 
responding canonical and phantom fields possess O(A^) symmetries. Assuming exponential 
potentials wc find that this 0{N) quintom paradigm exhibits novel properties comparing 
to the simple canonical and phantom scenarios. In particular, we find that the universe 
cannot result in a quintessence- type solution with w > —1, even in the cases where the 
phantom field seems to be irrelevant. On the contrary, there are always late-time attractors 
which correspond to accelerating universes with w < ~1 and with a recent crossing of the 
phantom divide, and for a very large area of the parameter space they are the only ones. 
This is in contrast with the previous simple-quintom results, where an accelerating universe 
is a possible late-time stable solution but it is not guaranteed. 
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1 Introduction 



Many cosmological observations, such as SNe la [jj, WMAP [2], SDSS [3], Chandra X- 
ray observatory [3], etc., reveal that our universe is undergoing an accelerating expansion. 
In addition, they suggest that it is spatially flat, and consists of about 70% of dark energy 
with negative pressure, of 30% of dust matter (cold dark matter plus baryons), and of an 
negligible amount of radiation. Although the nature and origin of dark energy could perhaps 
understood by a fundamental underlying theory unknown up to now, physicists can still propose 
some paradigms to describe it. The most obvious theoretical candidate for dark energy is the 
cosmological constant [SlEllTj which has the equation of state w = —1. However, it leads to the 
two known difficulties [8] , namely the "fine-tuning" problem (why is the current vacuum energy 
density so small), and the "cosmic coincidence" one (why are the densities of vacuum energy 
and dark matter nearly equal today since they scale very differently during the expansion 
history). 

There have been many efforts to resolve these problems [6l[9], but none could offer a robust 
and undoubted solution. This fact led many theoretical physicists to construct alternative 
frameworks, such is the dynamical dark energy scenario, by assuming that the vacuum energy 
is cancelled to exactly zero by some unknown mechanism, and introducing a dark energy 
component with a dynamically variable equation of state. The dynamical dark energy proposal 
is often realized by some scalar-field mechanism which suggests that the energy form with 
negative pressure is provided by a scalar field evolving downwards a proper potential. A 
large class of scalar-field dark energy models have been studied in the literature, including 
quintessence lLO\ , K-essence [H] , tachyon , phantom p^j , ghost condensate [H] , holographic 
dark energy [15], bulk holographic dark energy [16] and many others. 

A primary scalar field candidate for dark energy was the quintessence scenario |10j . which 
consists of a fluid with equation-of-state parameter lying in the range —1 < w < —1/3. On the 
other hand, for the phantom model [13] of dark energy, which consists of a scalar field with a 
negative sign of the kinetic term in the Lagrangian, one always obtains w < —1. Thus, neither 
the quintessence nor the phantom alone can fulfill the transition from t(;>— ltow<— 1 and 
vice versa, although the cosmological observations mildly favor models where such a transition 
was indeed realized and in particular with w crossing —1 in the near past. As it was indicated 
in the literature [17], the consideration of the combination of quintessence and phantom in a 
unified model, leads to the fulfillment of the aforementioned transition through the w = —1 
divide. This model, dubbed quintom, can produce a better fit to the observational data. 

The generalization of quintessence and phantom models to fields with O(A^) symmetry have 
been performed in [18] and [19] respectively. As it was shown, the behavior of the corresponding 
dynamical systems, in specific areas of the parameter space, can be qualitatively different than 
then single-field models. In the present work we are interested in investigating the generalized 
quintom model with an 0{N) symmetry. As a specific potential form we impose the exponential 
dependence on the corresponding fields, since exponential potentials are known to be significant 
in various cosmological models [20]. We perform a complete phase-space stability analysis of 
the corresponding autonomous system and we extract its attractor properties. We find that 
the O(A^) quintom has qualitatively novel properties compared to the corresponding simple 
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phantom, quintessence and quintom models. In particular, cosmological solutions with w < —1 
are the only attractors for a very large area of the parameter space, and the crossing through 
the w = —1 divide has a large probability to be realized. 

The plan of the work is as follows: In section[2]we construct the quintom model with O(A^) 
symmetries and in section [3] we perform its complete stability analysis. In section H] we discuss 
the cosmological implications of our results, and finally section [5] is devoted to conclusions. 



2 0{N) quintom 

We consider a flat Robertson- Walker metric: 

ds^ = dt^ -a^{t)diL^. (1) 
The Lagrangian density for a quintom model with O(A^) symmetries is: 

L = \g^^ [(a^<l>")(5,cl>°) - {dy){d,a^)\ - V^{\^^\) - V.{\A), (2) 

where is the component of the canonical field, with a = 1, 2, • • • , A'^^, and is the compo- 
nent of the phantom field, with /? = 1, 2, • • • , N^. Note that in the general case, the dimension- 
ality of the multiplets of the two fields, is not the same. Fortunately, these dimensionalities do 
not appear in the final form of the equations. 

In order to impose the O(A^) symmetries, following [18], we write: 

<I>^ = COS99$i(t) 

= i?$ {t) sin 99$^ (t) COS (/?$2 (i) 

= sin (t) sin (/?$2 (*) cos (/3$3(i) (3) 



and similarly: 



^ = i?$(t)sin(/j$j(t)---sin(/j$^_^_2(0cos(/7$j^^^_-^(t) 
= i?$(i)sin(/9$j(t)---sin(/9$ 2(t)sin99$^ (t), 



cr^ = Ra{t) cos Lp a A^) 

= Ra (t) sin ip^^ (t) cos ip^z (i ) 
cr^ = Rait) sin ip^^{t) sin ip^^{t) cos ipa.^{t) (4) 



cr^"^ ^ = Ra{t)sinipa^{t)---sinipa^^_^{t)cos(paj^^_^{t) 
a^" = Ra{t)sinipa,{t)---sin<fa^^^2{t)sin<fa^^^^{t). 

Thus, we have explicitly used the properties = R^ and \a^\ = R^- Furthermore, we 
assume that the potentials V$(|$"|) and Vo-df^l) depend only on and R^j respectively. 
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The action for the universe is as usual 

S 



where g is the determinant of the metric tensor g^^, Rg is the Ricci scalar and G is the Newton's 
constant (in the following we will instead use = SttG) . is the pressure of the barotropic 
fluid which constitutes the matter content of the universe, with equation of state = {'~f—l)p^, 
with the constant 7 in the interval < 7 < 2. 

The Einstein equations for the angles ip^^ and can be easily derived through the 
corresponding variation of the action. However, they are irrelevant for the purpose of this 
work since we are going to use only the radial equations for the fields, i.e. the equations 
determining the evolution of R^{t) and Ra{t), plus the Friedmann equations. Hence, we have: 

^2 

-ff^ = Y [p-y + /?$ + Pa] , (6) 



a\ / 37 



a J 3 



■ (7) 



2 Ij p^ + 2p<s, + 2pa + V^{R^) + Va{Ro 
H,,3™,-i|-^!MM.O. (8) 

In equations ([6|),([7]), the energy density and pressure of the canonical and the phantom fields, 
are given by: 

P'i> = l{Rl + ^) + V^{R^) 



p^ = l^Rl + ^)-V^{R^) (10) 



and 



P^ = -l{Rl + ^)+Va{Ra) 

Pa = -l{Rl + ^)-Va(,Ra). (11) 

In addition, the effect of the "angular component" of the system is embedded in the radial 
equations dS]),® as an effective term containing the constants $1$ and 0,^, which are determined 
by the values of the first integrals of motion [18]. Finally, H is Hubble parameter. 
The equation of state for the O(A^) quintom is: 

+ „ Rl-Rl + ^(§-§)- 2mR^) + VaiRa)] 
w = = ) " . (12) 



Rl-Rl + ^[-^-Wj+ 2[^$(i?$) + VARa)] 
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Thus, the constructed quintom model could produce a value w < —1 ii 

1 fnl nl 



Rl + 



^1 



<2\V^{R^) + V„{R„] 



(13) 



3 Stability analysis of the O(A^) quintom 

We are interested in investigating the attractor properties of the O(A^) quintom model, impos- 
ing exponential potentials, since they are known to be relevant in various cosmological models 
j. In particular we consider: 



V$(i?$) = V$(j exp(— AKi?$) 
Va{Ra) = exp(-AKi2^). 

Using this specific potential ansatz, the radial equations of motion 

+ 3HR^ - -g^g- - \k V^{R^) = 0, 



become: 



Rfj + 3HRq 



a^Ri 



\kV„{R,) = Q. 



(14) 



(15) 



(16) 



Furthermore, using the definitions for the energy densities and pressures (|1U|) . (|11|) . the Fried- 
mann equations ([I]),© can be re- written as: 



H 



K 

y 



Ri 



a^Rl 



+ V^{R^) 



a^Rl 



+ V^{R,) 



(17) 



(18) 



Finally, the equations close by considering the evolution of the barotropic (matter) density: 

p^ = -3H{p^+p^). (19) 

In order to perform the stability analysis of the 0{N) quintom model, we have to trans- 
form the dynamical system ()15p - (ll8p into an autonomous form ^21j. This will be achieved by 
introducing the auxiliary variables: 



V6H 



R^ 



K^V^jR^) 

^/6H a3R<i> 
kR^ 



V6H 



Ra 



K^VaiR.) 

1 



KRfj 



(20) 
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together with M = log a. 

Using these variables, we result in the following autonomous system: 

dM 

dXa 

dM 

dm 

dM 

dM 
dz^ 

dM 

dZa_ 

dM 









3xa + VqzI^o- - \ \>^yl 

V z 


\v^T, - 




i^y.T, - 




Iz^T, - 


32$ — ^/Qx^z^^i^ 


^z^n - 


3Za - VeXaZaCa 



dM 



-^/6e>., (21) 



where Ti = 7 (1 - x| - y| - z| + - _|_ _^ 2(x| + z| - - z^). Note also that the 
Friedmann equation ()18p leads to the constraint equation: 

2 2 I 2 I 2 I 2 2 I ^^^7 . 

- + y$ + + - + = 1. (22) 
Finally, in terms of the auxiliary variables, the equation of state for the quintom ()12p becomes: 

^y>2 ^2 , 2 ^ 2 _i_ ^2 ^2 

^ "''ja '^Q- Ha ~r ^cr (23) 

The critical points {x^c,Xac,y'^c,ycTc,z^c,Zac,S,^c,S,cTc) of the autonomous system ([2T]1 are 
obtained by setting the left hand sides of the equations to zero. The real and physically 
meaningful of them are presented in table [H 

In order to determine the stability properties of these critical points, we proceed as follows: 
We expand the auxiliary variables around the critical points as 













= y<j.c + v<!, 


ya 


= yac + Va 




= Z^c + 


Za 


= Zac + Wa 




= C*c + X* 




= C<7C + Xa- 



(24) 
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Table 1: The real and physically meaningful critical points of the autonomous system (I2ip . 

In expressions ([MI), Ui^,Ua, v^,Va,w^,Wa, Xo- are just the perturbations of the variables near 
the critical points and we consider them forming a column vector denoted as U. Inserting these 
expansions into the autonomous system ()21|) . we can obtain the equations for the perturbations 
up to first order as: 

U' = Q • U, (25) 

where the prime denotes differentiation with respect to M. 

For the critical points {x^c,Xac,y^c,yac, z^c, Zo-c,(,^c,Cac), the coefficients of the perturba- 
tion equations form a 8 x 8 matrix Q. Thus, for each critical point of table [H the eigenvalues 
of Q determine the type and stability of this specific critical point. The explicit form of the 
matrix Q and its eigenvalues for the critical points of table [H i.e for the autonomous system 
(j21|) . are given in the Appendix. In table [2] we present the results of the stability analysis. 
In addition, for each critical point we calculate the values of w (given by relation ()23|] ). and 
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of (given by the constraint (j22p ). Thus, = means that the universe is dominated 
completely by the quintom fields. In the next section we discuss the cosmological implications 
of the obtained results. 

4 Cosmological implications 

The critical points K and L are stable nodes for every 7 and every A and thus they constitute 

2 

late-time attractors. Since = they correspond to quintom-dominated universes, with 
the quintom equation-of-state parameter being w = —1 — The fact that x^c = y$c = 



7 





r lYl pn OP 
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A 


V 7 and V A 


Stable node for \x„ > 
Unstable node for -y^e - \^ < \x„ < y^a- 
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-1 
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u 


T3 
D 


— w — 

V 7 and V A 
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1 
i 


U 




U < A < 


Unstable node 


1 
i 


U 


L) 


U < A < 
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U 


"P 

rj 
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Saddle point 


1 _L \^ /Q 

— i + A /O 


n 
U 


"P 
r 


A < D 


Saddle point 


1 _l_ \2 /q 
— i + A /O 


U 




V 7 and A ^ 


Saddle point 


7 - i 


i - ^7/ A 


H 


V 7 and A 7^ 


Saddle point 


7-1 


1-37/A^ 


I 


V 7 and A 7^ 


Saddle point 
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V 7 and A 7^ 
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V 7 and V A 
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-I-A73 
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V 7 and V A 


Stable node 


-1-AV3 






Table 2: The properties of the critical points of the autonomous system ([2T]) . 



implies that the phantom component of the quintom plays the main role in this solution. Thus, 
this solution corresponds to the phantom solution of the literature, with the expected u; < — 1. 
However, the more complex dynamics of the quintom model makes the system qualitative 
different from the simple phantom one, even in this particular solution where the canonical 
field seems to have a trivial contribution. Indeed, in the simple phantom case [19] a stable 
node exists only when A^ < 6, which in turn leads to a lower bound for w. As we see this is not 
true in the 0(A'^) quintom paradigm, where the stable node is independent of the parameter 
values. Thus, although in the previous studies an accelerating universe was not always a late- 
time attractor, in the present analysis we show that such a cosmological era is always a stable 
equilibrium solution. In addition, the possible values of w are not bounded. 

The critical points A and B are a novel property of the present quintom model. First of all, 
since for every there is always a x$ which lead to zero left hand sides of the autonomous 
system (12 ip (given the nullification of the other variables), A and B corresponds to loci of 
critical points. They correspond to quintom-dominated solutions, with w = 1. Note that in 
this case both the canonical and the phantom field are relevant for the cosmological evolution. 
The critical point B is unstable. However, in the case Xxa > \/6, the critical point A is a stable 
node. This behavior is different form both the simple canonical [20] and the simple phantom 
[19] fields. 

The critical points C, D and I correspond to quintom-dominated solutions, and especially 
to solutions where the constraint ()22p is dominated only by the kinetic energy of the canonical 
and the phantom fields. As expected, these "kinetic" -dominated solutions are unstable and 
are only expected to be relevant at early times. Furthermore, the critical point J corresponds 
to a quintom-dominated unstable solution. 

The critical points G and H correspond to saddle points of the dynamical system at hand. 
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for every 7 and every A 7^ 0. Note that although for a specific combination of these parameters 
(7 > I and > -^p^) two of the eigenvalues form a complex conjugate pair with negative 
real part, the corresponding critical points are not stable spirals but saddle points, due to 
the presence of a positive eigenvalue. This is a crucial difference comparing to the simple 
canonical field [20j and reveals that the more complex quintom dynamics changes qualitatively 
the behavior of the system, even if the phantom field does not seem to play an important role 
in this case {x^c = Uac = 0). The value of w for the critical points G and H is 7 — 1, thus for 
the allowed 7 values it corresponds to —1 < u; < 1. In other words, this solution corresponds 
to the quintessence universe and the fact that is not a stable point, i.e. a late-time attractor, 
means that in the model at hand the universe cannot result in a quintessence-type solution 
with w > —1, even in the cases where the phantom field seems to be irrelevant. Finally, we 
mention that for these points = 1 — Thus, although G and H exist for every 7 and 
every A 7^ 0, the corresponding solution has a physical meaning only when A^ > 87. This fact 
was also mentioned in [20] for the simple canonical model. 

The critical points E and F, are saddle points. They correspond to quintom-dominated 
universes, with w = — 1 + The fact that x^-c = Vac = means that the phantom field is not 
relevant for this solution, and the cosmological evolution is driven by the canonical field. Thus, 
this solution also corresponds to the quintessence universe. However, while this solution is a 
stable node, i.e. a late-time attractor, for the simple quintessence model in the case A^ < 87 
|20j . in the present quintom model it is always a saddle point. Therefore, the quintessence- 
type solution cannot be an equilibrium point for the cosmological evolution of the universe. 
Once again we see that the more complex 0{N) quintom dynamics implies qualitatively new 
behavior for the cosmological evolution comparing to the previous simple models. 

From the aforementioned investigation we conclude that the 0(iV) quintom has a variety of 
critical points. It is interesting to see that for not very steep potentials and/or a small x^, the 
attractor A disappears. Thus, the attractors K and L are the only ones for a very large area 
of the two-dimensional parameter space (7, A). And since the initial value of w can be easily 
above —1 we conclude that the transition through the w = —1 divide has a large probability 
to be realized, as the system is attracted be the attractors K and L. 

5 Conclusions 

In the present work we investigated the quintom paradigm of dark energy in the generalized 
case where the corresponding canonical and phantom fields possess 0{N) symmetries. Our 
analysis reveals that this O(A^) quintom model exhibits novel properties comparing to the 
simple canonical and phantom scenarios. In particular, we find that the universe cannot result 
in a quintessence- type solution with w > —1, even in the cases where the phantom field seems 
to be irrelevant. On the contrary, there are always late-time attracting equilibrium points 
which correspond to quintom-dominated, accelerating universes with w < —1, and for a very 
large area of the parameter space they are the only ones. This is in contrast with the previous 
simple-quintom results, where an accelerating universe is a possible late-time stable solution 
but it is not guaranteed. In addition, the more complex quintom dynamics, refutes the lower 
bound of w obtained previously in the simple phantom paradigm, making its possible values 
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unlimited. Finally, the transition form above to below of the phantom divide w = —1, has 
a large probability to be realized. These features make the present 0{N) quintom model 
consistent with cosmological observations. 



Appendix: stability of the critical points 

For the critical points {x^c,Xac,y^c,yac,Z'^c,Zac,S,^c,S,cTc), the coefficients of the perturbation 
equations form a 8 x 8 matrix Q, which for the autonomous system ()21|) reads: 



Qii 

Ql2 
Ql3 
Ql4 

Q21 

Q22 

Q23 
Q24 

Q31 

Q32 

Q33 
Q34 
Q41 

Q42 

Q43 

Q44 

Q55 
Qee 



-1 + Y+4c(2-7) 
y$c(\/6A - 37x$c) 

-37X$c Vac 

2-7) 



3 X^c Xcrc 



-1 + Y+4c(7-2) 



Uaci-V^X - 37X0 



-A + 3x$c(2 

3 Xac (7 - 2) 
-^-Xx^, + — - 

-37y*c Vac 

3 X$c Vac (2 - 7) 

3" 

A + 3 x^c (7 
3T2 



7) 



2 



Q 2 

3 7y$c 



2) 



o 2 



-1 + 



T2 



^2 

-1 + — 

2 



where T2 = 7(1 



r-2 

"(TC 



vie - vie) + 2(a:L 



All the other components of Q are 



zero for the physically interested critical points of tabled! 

The eigenvalues z/j (i = 1, • • • , 8) of Q for the specific critical points of tabled! i.e for the au- 
tonomous system ([2T|) , are presented in tabled where T3 = ^^A2(2 - 7)(2472 + 2A2 - 97A2). 
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Thus, by determining the sign of the real parts of these eigenvalues, we can classify the corre- 
sponding critical point |21j . 
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1(7-2) 
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!(7-2-T3) 


i(7-2 + T3) 


1(7-2) 


1(7-2) 


1(7-2) 
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6-37 


3- ^|V6 + A2 
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6-37 


3+ y|V6 + A2 
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- 3 
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2 
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A^ 
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Table 3: The eigenvalues of the matrix Q of the perturbation equations of the autonomous 
system (f2T]) . 
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